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ABSTRACT 

The microquasar GX 339—4, known to exhibit powerful compact jets that dominate its radio 
to near-infrared emission, entered an outburst in 2010 for the fifth time in about fifteen years. 
An extensive radio to X-ray multi-wavelength campaign was immediately triggered, and we 
report here on ESO/FORS2+ISAAC optical and near-infrared spectroscopic observations, 
supported by ATCA radio and RXTEISwift X-ray quasi-simultaneous data. GX 339—4 was 
observed at three different epochs, once in the soft state and twice in the hard state. In the soft 
state, the optical and near-infrared continuum is largely consistent with the Raleigh- Jeans tail 
of a thermal process. As an explanation, we favour irradiation of the outer accretion disc by its 
inner regions, enhanced by disc warping. An excess is also present at low frequencies, likely 
due to a M subgiant companion star. During the first hard state, the optical/near-infrared 
continuum is well-described by the optically thin synchrotron emission of the compact jet 
combined with disc irradiation and perhaps another component peaking in the ultraviolet. 
The spectral break where the jet transits from the optically thick to thin regimes, located 
below 1.20 x 10 14 Hz, is not detected and the extension of the optically thin synchrotron 
is consistent with the 3-50 keV spectrum. In contrast, the emission during the second hard 
state is more difficult to understand and points toward a more complex jet continuum. In 
both cases, the near-infrared continuum is found to be variable at timescales at least as short 
as 20 s, although these variabilities are smoothed out beyond a few hundred seconds. This 
implies rapid variations - in flux and frequency - of the location of the spectral break, i.e. 
dramatic short timescale changes of the physical conditions at the base of the jet, such as the 
magnetic field and/or the base radius. 
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1 INTRODUCTION 

While microquasars radiate over the whole electromagnetic spec- 
trum, the optical and near-infrared (near-IR) domains are of strong 
interest. Indeed, most of the components in the system can be 
detected at these wavelengths, including the companion star, the 
accretion disc, and material ejecta. In the last category, the best 
known are the so-called compact jets, which are only detected 

* Based on observations performed with European Southern Observatory 
(ESO) Telescopes at the Paranal Observatory under programmes ID 284. D- 
5056 and 285.D-5007 
t E-mail: frahoui@cfa.harvard.edu 



in the hard state (HS) and emit both through optically thick 
and thin synchrotron radia t ion. Within the most acce pted scheme 
dBlandford & Konigll979l ; lFalcke & BiermanrJI 19951) , their emis- 
sion is well-modelled by a flat or weakly-inverted power law (op- 
tically thick regime, F„ oc v a with —0.1 ^ a ^ 0.7) from 
the radio to some spectral break, beyond which a ranges between 
—0.4 and —1 depending on the electron energy distribution (op- 
tically thin regime). The location of the break, thought to occur 
in the IR domain, is a crucial piece of information as it is closely 
related to the physical conditions at the base of the jet, such as 
the magnetic field, the base radius of the jet, and the total en- 
ergy of the electron population. So far, it has only been detected 
in three sources, GX 339-4, 4U 0614+091, and Cygnus X-l 
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Table 1. Summary of all the observations of GX 339—4 we made use of in this study. We give the observation number, the day of observation with FORS2 
and ISAAC (in MJD), the day of quasi-simultaneous coverage with RXTE, Swift, and the ATCA telescope, as well as the ATCA flux levels at 5.5 and 9.0 GHz 
(in mJy). 



Obs. # 


FORS2 


ISAAC 


RXTE 


Swift 


ATCA 


ATCA flux density 
5.5 GHz 9.0 GHz 


1 


55261.39 


55261.37 


55261.19 


55261.05 


55261.89 


9.02±0.10 


9.60±0.05 


2 


55262.39 


55262.38 


55262.93 




55262.91 


8.24±0.05 


8.05±0.10 


3 


55307.29 


55307.24 


55307.02 











dCorbel & Fended |2002|; IGandhi et alJliout iMigliari et aUl2006l. 
l2010l :l Rahoui et alj|2011f) . with additional obs ervational evidence 
that this break can move with lum i nosity dCoriat et alj 1 20091 : 
iRahouietal .1120 111 ; IGandhi et al]201 ll ; iRussell et alj201 lh . 

If a broken power law is a good approximation of a compact 
jet, this very simple modelling is also plagued with severe limita- 
tions. First, it is doubtful that the jet can sustain a single power law- 
like spectrum from the radio to the spectral break, and a change of 
the optically thick spectral index is expected. Moreover, their emis- 
sion could simply be more complicated than a br oken power law, 
in pa rticular when cooling is taken into account dPe'er & Casellal 
[20091) . Finally, compact jets could contribute to the optical/near- 
IR (OIR) through other processes such as pre-shock synchrotron, 
which would en hance significantly th eir emission at these wave- 
lengths (see e.g. iMarkoff et alj |2005). OIR spectroscopic obser- 
vations of microquasars are therefore crucial for correctly con- 
straining the properties of compact jets; GX 339—4 is obviously 
one of the best targets to fulfil this goal a s the jets strongly domi- 
nate the OIR emiss ion in the HS (see e.g. lBuxton & B ail vrj| 20041 ; 
iHoman et ai]|2005h . 

In early 2010, GX 339-4 entered an outburst, only the fifth 
time in fifteen years. Considering the rarity of such an event, we 
triggered a large multi-wavelength observational campaign of the 
source, which included X-ray (INTEGRAL, RXTE, Swift), OIR 
(ESO, Faulkes), and radio (ATCA). Here, we report on the OIR 
spectroscopy of the source in the HS and soft state (SS), focusing 
on the spectral continuum in the presence and absen ce of com- 
pact jets. A companion paper dCadolle Bel et alj|201 lh presents an 
overview of this multi-wavelength campaign, with a particular fo- 
cus on X-ray behaviour. Finally, a detailed study of the spectro- 
scopic content of the OIR spectra will be object of a third paper 
(Rahoui et al., in prep). Sect. 2 is devoted to the data reduction and 
analysis. In Sect. 3, we present spectral energy distribution (SED) 
modelling. We discuss the outcomes in Sect. 4 and conclude in 
Sect. 5. 



2 OBSERVATIONS AND DATA REDUCTION 

GX 339—4 was observed through low-resolution spectroscopy on 
2010 March 6 (hereafter Obs. 1), March 7 (Obs. 2), and April 21 
(Obs. 3) quasi-simultaneously in the optical with ESO/FORS2 and 
in the near-IR with ESO/ISAAC, mounted on the VLT/UT2 and 
UT3, respectively (PI F. Rahoui). In the framework of the Galac- 
tic bulge monitoring, the source was also quasi-simultaneously 
observed with the Rossi X-ray Timing Experiment (RXTE, the 
data are immediately public). Moreover, the Swift satellite also 
observed the field of GX 339—4 on a regular basis during 
the 2010 outburst, and one of these pointed observations (Obs. 
ID 00030943010, P.I. M. Cadolle Bel) was conducted quasi- 
simultaneously to Obs. 1. Finally, GX 339—4 was detected in the 
radio domain during Obs. 1 and 2 with the Australia Telescope 
Compact Array (ATCA, P.I. S. Corbel) at 5.5 and 9.0 GHz. All 
these observations are listed in TableQ] 



2.1 Near-IR ESO/ISAAC observations 

GX 339—4 was observed through the SW— LR1 arm (l"slit aper- 
ture) with the J, H, and K filters (R ~ 600 and 1.10-2.5 fim 
spectral coverage). The sky was clear, and the airmass and see- 
ing at corresponding central wavelengths were lower than 1.2 and 
OV75, respectively. The total exposure time in each filter was 240 s, 
divided in 4 x 3 x 20 s dithered frames for standard ABBA sky 
removal. The same nights, the telluric standard stars HIP 84086 
(B8/9V, March 6 2010), HIP 86735 (B5/6III, March 7 2010), and 
HIP 68124 (B2/3V, April 21 2010), were observed under the same 
conditions. 

We reduced all the data with IRAfQby performing correction 
for crosstalk and bad pixels, flatfielding, and sky subtraction. The 
spectra were then extracted, wavelength-calibrated by matching 
the spectral dispersion axis with that of an ArXe arc, and the atmo- 
spheric signatures were corrected with the telluric task. This 
routine aligns the source's and standard star's spectra and scale 
them according to their respective airmass; they are then divided to 
remove the telluric features and atmospheric extinction. The out- 
put, once corrected for the exposure time discrepancies between 
the source and the standard star observations, can therefore be mul- 
tiplied by the flux-calibrated standard star's spectrum to obtain that 
of the source. Unfortunately, very few spectro-photometric stan- 
dard stars exist in the near-IR, and none of them was observed 
during our program. The way we proceeded to flux-calibrate the 
GX 339— 4's ISAAC spectra therefore deserves some detailed ex- 
planations. 

We started from three synthetic Kurucz spectra corresponding 
to the HIP 84086, HIP 86735, and HIP 68124 spectral types. They 
were then reddened, convolved by Cousins-Bessel and 2MASS fil- 
ters, and scaled to match, through \ 2 fitting, the measured B, V, J, 
H, and K flux densities of our three standard stars (see Figure lATl . 
The initial value of the optical extinction Av along their line-of- 
sight (LOS) were obtained from their colour indices using the stan- 
dard relation Av — Rv x E(B — V) with Rv = 3.1 (correspond- 
ing to the diffuse ISM) and E(B — V) = (B — V ) - (B - V) 
The intrinsic colours (B — V)o were taken from ICramer] dl984l) 
for the corresponding spectral types, an d we made use of the ex- 
tinction law given in iFitzpatrickl d 19991) to get the A\ values at 
other wavelengths. Once flux-calibrated, the standard stars' syn- 
thetic spectra were multiplied by the telluric outputs to get 
the final GX 339— 4's flux-calibrated spectra. It is worth men- 
tioning that the ESO/REM J, H , and K photometric ob servations 
quasi-simultaneous with Obs. 1 dCadolleBeletalj201lh . are con- 
sistent, with less than 5% discrepancies, with our result. Therefore, 
in our analysis, we adopt an intrinsic 5% uncertainty for the flux- 
calibration. 



1 I RAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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Table 2. Best parameters obtained from the fits to the GX 339—4 X-ray spectra during Obs. 1, 2, and 3, with the model 
PHABS X REFLECT X (DISKBB+GAUS SI AN+POWERLAW). The errorbars are given at the 90% confidence level. 



Obs. # 


7V H 


Q/2tt 
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r 




X? (dof) 




(X10 22 cm" 2 ) 




(keV) 


( V°> ^ 
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(keV) 


(keV) 












Vc/lflkpc) 








1 




4Q+0.10 

u -^ a -0.09 


n 9 ,+0.03 
u - zo -0.02 


3.50^ 


X 10 4 


6 31+° 19 
°-° -0.23 


fl fil+0-35 
U - 8i -0.30 


-, 7 o+0.02 
- L - lo -0.02 




1.06 (651) 


2 


0.58 (fixed) 


39+ 11 








6 24+ - 18 
D - z *-0.22 


62+ 34 


i 79 +0.03 
1 - ' -0.03 




0.95 (74) 


3 


0.58 (fixed) 




°.83ig:E 


o rio + 0.88 


X 10 3 


5.97+0.44 


1 qq+0.48 
1 -°°-0.31 


Z - 4D -0.16 




1.03 (73) 



2.2 Optical ESO/FORS2 observations 

We used the 300V and 3001 grisms (1" aperture) with the GG435 
and OG590 filters, respectively. The exposure time was set to 120 s 
in each filter, and the atmospheric conditions were good, with a 
clear sky and airmass and seeing lower than 1.14 and Of.'l, respec- 
tively. For flux calibration purposes, the ESO spectro-photometric 
standard stars Hiltner 600 (B3V) and CD-32 9927 (A0V) were 
observed under the same conditions. We reduced the data using 
the I RAF suite. The process consisted of bad pixels correction, 
bias subtraction, flatfielding, spectral extraction, and wavelength 
calibration by matching the spectral dispersion axis with that of 
an HeAr arc lamp. The resulting GX 339—4 ID spectra were then 
flux-calibrated using the CTIO extinction and flux tables0 avail- 
able with I RAF. The flux levels in the overlapping region of the 
GG435 and OG590 filters appear consistent, with discrepancies 
less than 5%, which allowed us to safely combine them. How- 
ever, the strong limitation of this method is that the resulting op- 
tical spectra do not overlap with the near-IR ones because the 
ESO spectro-photometric standard stars are not defined beyond 
1.02 /xm. To cover the gap between the optical and near-IR, we 
also reduced and flux-calibrated the FORS2 spectra with the same 
procedure as that used for the ISAAC data (see Figure IA2I for 
the flux-calibrated synthetic Kurucz spectra). This gives very sim- 
ilar results, as showed in Figure [BJ] (left panel). The optical spec- 
tra flux-calibrated with the ESO spectro-photometric standards are 
displayed in green, those flux-calibrated with synthetic spectra in 
red. We also display the J band spectra in blue, and a good match 
between the optical and near-IR continuum levels is reached (right 
panel). Note that the displayed red spectra were rebinned to a third 
of their resolution to improve the visibility of the optical/7 overlap- 
ping region. However, only the original ones were used for SED 
fitting. 



2.3 RXTE and Swift observations 

We present here the reduction and analysis of the high-energy data 
we made use of. We use this only for the determination of the spec- 
tral state of GX 339-4 during the FORS2+ISAAC observations, 
as well as the level of emission from the main components, i.e. 
the accretion disc and the corona/jet. A far more comprehensive 
X-ray study of GX 339— 4 during the 2010 outburst is presented in 
ICadolle Bel et alj J201 lh . 



2.3.1 Data reduction 

The RXTE data were reduced using the HEASOFT v. 6.11 software 
package, following the standard procedure described in the RXTE 
cookbook Spectra were extrac ted from the Proportional Counter 
Array (PCA: |jahoda et alj|2006l) . in the 3-50 keV range. We only 



2 ctioextind.dat, h600.dat, wAcdi2.dat 

3 http://heasarc.gsfc.nasa.gov/docs/xte/data_analysis.html 



used the top layer of the Proportional Counter Unit (PCU) 2 as it is 
the best-calibrated detector out of the five PCUs. Systematic errors 
of 0.5 per cent were added to all channels. 

We reduced the data from the X-ray telescope (XRT) on board 
Swift with the xrtpipeline v. 0.12.6 using standard quality 
cuts and event grades 0D2. The data were collected in windowed 
timing (WT) mode to avoid pile-up. We used xselect v. 2.4 to 
extract source and background spectra using square boxes as ex- 
traction region, in the range 0.5-10 keV. The boxes have a size of 
40 pixels along the only spatial dimension available in WT mode. 
We generated the ancillary response file (ARF) with xrtmkarf 
and used the latest version (v. 012) of the response matrices pro- 
vided by the Swift team. We finally rebinned the spectrum to get a 
minimum of 20 counts per channel. 



2.3.2 Spectral state classification 

To identify the spectral state of GX 339—4 during Obs. 1, 2, and 
3 we fitted, with XSPEC v. 12.7, the three high-energy spectra 
with a model combining the con tributions from the accretion 
disc (DISKBB. [Mitsuda et al]|l984l) . an iron line around 6.4 keV 
(GAUSSIAN), as well as a power law, all modified by photo- 
electric extinction ( P HABS, a bundances and cross-sections from 
Anders & Grevessel 1 19891 ; iBalucinska-Church & McCammonl 
19921) . When needed, a reflection component REFLECT was added 
for a n inclination angle fi xed to the most recent measured value, 
46° dShidatsuetal]|201ll) . However, it is worth remembering that 
se veral other measurem ents between 10 and 50° exist (see Table 6 
in lShidatsu et al.|[2.01 lL for a review), but the use of a 10° incli- 
nation angle does not change significantly the fits. Moreover, the 
disc's contribution during Obs. 2 was too low to be detected with 
PCA so that DISKBB was not used in the modelling. The best-fit 
parameters are listed in Table|2] The results unambiguously show 
that GX 339-4 is in the HS during Obs. 1 and Obs. 2 (hereafter 
HS1 and HS2) , and in the SS during Obs. 3 (hereafter SSI). 



3 NATURE AND PROPERTIES OF THE OIR EMISSION 

We use spectral fitting to understand the origin of the OIR con- 
tinuum of GX 339—4. To that purpose, we built the SEDs for the 
three observations with all the available data listed in TableQ] The 
OIR spectra were c orrected for ISM e xtinction along the LOS with 
the law derived in iFitzpatrickl jl999h . Among the several consis- 
tent values for the optic al absorption (see e.g iBuxton & Vennesl 



120041) . we adopted Ay = 3.25 ± 0.50 (see 



20031: iHvnes et al 

Gandhi et alj|201 1 , and reference therein), but it is however worth 
remembering tha t it is not well constrai ned and may be higher (e.g. 
Ay = 3.7±0.3, lzdziarski et al.ll998l) . In adopting the lower Ay. 
our subsequent results, especially those related to the disc emission 
in the optical, should be considered as lower limits. Throughout 
HS1 and HS2, the source was detected in the radio domain with 
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Figure 1. Near-IR sub-spectra for HS1 (top), HS2 (middle) and SSI (bottom), and in each band (K to J, left to right) of GX 339—4. Each sub-spectrum 
corresponds to a 20 s integration. 



a relatively flat emission and we therefore expec t a strong contri- 
bution from the compact jets dCorbel et alj|2003l) . In contrast, the 
continuum should be dominated by the accretion disc during SSI. 



3.1 Near-IR continuum variations at short timescales 

Figure [T] displays, for the three observations and in each near-IR 
band, the twelve deredenned "sub-spectra" that were all obtained 
with a good signal-to-noise ratio (S/N). The comparison of their re- 
spective continuum may therefore give reliable information about 
the intrinsic spectral variations of GX 339—4 at timescales at least 
as short as 20 s, both in the HS and SS. However, the optical spec- 
tra being the result of one exposure only, their time-dependent be- 
haviour cannot be known. 

During SSI, the near-IR emission of GX 339—4 is constant 
in each band over the 12 x 20 s exposures. On the contrary, in 
the HS the emission is strongly variable in all bands on timescales 
as low as 20 s. While this behaviour is consistent with the near- 
IR fli ckering already observed through photometry l lCasella et al.l 
it could also be the result of changing observing conditions, 
in particular variable slit-losses due to a variable seeing. How- 
ever, the seeing was always lower than the 1" slit width, and the 



flux-calibration we reached duri ng Obs. 1 is consistent within 5% 
with the ESO/REM photometry dCadolle Bel et al.l201 ]]); this very 
likely points towards low slit-losses. Moreover, the atmospheric 
conditions were each time stable over the whole integration. In- 
deed, the airmass-corrected seeing at central wavelength, measured 
by ESO, was clustered in the ranges 0'.'62-0'.'73, 0'.' 63-0'.' 74, and 
0'.'57-0'.'66 during Obs. 1, Obs. 2, and Obs. 3, respectively. When 
considered per filter, the maximum-to-minimum variations were 
always lower than 10%, and even 5% during Obs. 1. This by itself 
does not exclude slit-loss variations, but it is reasonable to assume 
that their effect is weak. Finally, the fact that the continuum is con- 
stant during SSI and variable during HS1 and HS2 despite similar 
atmospheric and instrumental conditions strengthens the intrinsic 
origin of the variations in the HS. 

That said, the average maximum-to-minimum flux ratios 
/ma X //min in J, H, and K are 1.36 ± 0.08, 1.42 ± 0.08, and 
1.72 ± 0.13 during HS1, as well as 1.84 ± 0.12, 1.60 ± 0.10, 
and 2.29 ± 0.17 during HS2, respectively. Because J precedes H 
which itself precedes K, it is not relevant to compare the variations 
in a filter with respect to one another; it is however still possible 
to compare both hard states. In J and H, the continuum is clearly 
more variable during HS2. In K, although HS2 exhibits a higher 
maximum-to-minimum flux ratio, it only stems from three sub- 
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Figure 2. ISAAC+FORS2+PCA SED of GX 339-4 during SSI fitted with 
the irradiated accretion disc DISKIR model. Besides a relativistic iron line 
needed to fit an X-ray excess around 6 keV, a near-IR excess that we fitted 
with a ~3348 K black body is present at lower frequencies. We also display 
the irradiated disc with no OIR reprocessing. 



spectra and / ma x//mm = 1.15±0.04 otherwise. Overall, the vari- 
ations are therefore bluer in HS2 than in HS1, in agreement with 
their respective average continuum (see Sect. 3.3 and Figure [5); 
this likely points towards two different origins. 



3.2 The average OIR continuum during SSI 

The standard multicolour accretion disc alone is unable to re- 
produce the OIR spectrum, and a possible explanation is an ex- 
cess due to X-ray reprocessing within the disc. In order to get 
a quantitative unde rstanding of this phenomeno n, we considered 
the DISKIR model dGierlinskTeta l. 2008, 2009) combined with a 
Gaussian, to fit the OIR to X-ray SED. Roughly, it is an extension 
of the DISKBB model that includes comptonisation of the black- 
body disc's photons by a co rona of hot electrons (based on NTH- 
COMP. Izdziarski et al.ll996f) . as well as disc irradiation, both from 
the inner region of the disc itself and from the corona. The model 
has nine parameters: the disc's temperature fcT<ji SC and normalisa- 
tion N (same as DISKBB), the photon index V of the hard X-rays 
power law and the temperature kT c of the comptonising photons, 
the ratio between the corona's and the disc's luminosity L c /La, 
the fraction of hard X-ray emission that illuminates the disc fi n , 
the irradiated radius _f?i rr expressed in terms of the disc inner ra- 
dius, the fraction of soft X-ray emission which is thermalised in 
the outer disc / ou t, and the logarithm of the outer radius of the 
disc loglO(-Rout) expressed as a function of the inner radius. The 
first seven parameters are completely defined by the high-energy 
data while the two latter are characterised by the OIR ones. We 
fixed several parameters: (1) the comptonisation temperature kT c 
was unconstrained and pegged to the upper boundary whatever it 
was. We fixed it to 250 keV, a value high enough to mimic a power 
law like hard X-ray spectrum, (2) fi n was fi xed to 0.3, well-suited 
for the SS of microquasars ( 0.1 in th e HS. IPoutanen et al.lfl997l ; 
llbragimov et al.ll2005l : [Gilfanovll2010h (3) after several unsuccess- 
ful attempts, _ Rj rr was fixed to 1 . 1 as it systematically pegged to the 
lower value dGierlinski etal]|2009h . Moreover, this model alone 
was never able to describe the whole SED due to a near-IR excess. 
We therefore added a black body for completion, but it is important 
to note that this excess could have other explanations. The best-fit 
parameters are listed in Tableland the best-fit SED is displayed 



Table 3. Best parameters obtained from the 
fit to the GX 339-4 ISAAC+FORS2+PCA 
SED during SSI. The best-fit model is 
PHABSX(DISKIR+GAUSSIAN+BBODY). The er- 
rorbars are given at the 90% confidence level 



Parameters 



Best-fit value 



N H (xlO 



22 , 



-') 



fcT disc (keV) 

r 

kT c (keV) 

/in 

/out 

logl0(.Rout) (Rin) 

N" 
E ilon (keV) 

<T iron (keV) 

T BB (K) 
Rbb/Dbb (-R©/kpc) 

Xr (d-O.f) 



0.58 (fixed) 
0.78+™} 



+ 0.10 



2.38 u 
2.50 x 10 2 (fixed) 



19+ 003 
0.30 (fixed) 

0.0518$ 

4.301°:" 
4.35!°,;^ x 10 :i 
5.97 (fixed) 
1.33 (fixed) 
3.35t°-°* x 10 3 

1.181°;?? 
0.97 (4442) 



" Same as DISKBB 

in Figure [2] It is clear that accretion disc irradiation dominates the 
OIR continuum during SSI. Nevertheless, the derived proportion 
of reprocessed X-rays within the disc, / ou t = 0.051q'q2' ' s ver y 
high and at least one order of magnitude larger t han previously 
found in XTE J1817-33 dGierlinski et alj|2009h and twice the 
value for GRS 1915+105 dRahoui et al.ll2010h . Although we can- 
not exclude that this is due to an overestimation of th e optical ex- 
tinctio n, no lower value has ever been confirmed yet. Irlvnes et alj 
d2004l) derived A v ^ 2.64 from the interstellar Na doublet, but 
acknowledged that it was probably too low because the NaD fea- 
tures were saturated. It is therefore very likely that either another 
source of thermal emission contributes to the optical spectrum of 
GX 339—4 during SSI, and/or irradiation is enhanced. 



3.3 The average OIR continuum during HS1 

Figure[3]displays the deredenned OIR spectrum of GX 339—4 dur- 
ing HS1 (left panel), as well as the radio to X-ray SED (right 
panel). The OIR spectrum is best-described by two power laws 
whose spectral indices at a 3cr significance are am, = — 0.63 ± 
0.05, typical of optically thin synchrotron emission from compact 
jets, and av = 2.73 ± 0.37. Moreover the extension, from the ra- 
dio domain, of the optically thick emission (a ~ 0.15) matches the 
near-IR spectrum, in agreement with a spectra l break, which we d o 
not detect, located close to 1.2 x 10 14 Hz (see lGandhi et alj201lb . 
Likewise, the extension of the jet optically thin emission is consis- 
tent with the high energy spectrum beyond about 6 x 10 Hz, i.e. 
2.5 keV (see right panel, Figure[3]l. In the optical, the value of ay 
is larger than the spectral index expected for the Raleigh- Jeans tail 
of a thermal emission (~ 2), which may mean that besides disc ir- 
radiation, another component contributes to the ultraviolet/optical 
emission of GX 339—4 during HS1. Now, this could also be due 
to an overestimation of the optical extinction and considering the 
large optical error bars, a good fit is also obtained when av is fixed 
to 2 (Figure|4] left panel). In this case, am. = -0.74 ± 0.04, still 
consistent with optically thin synchrotron from the jet, but the ex- 
tension of the jet power law does not match the X-ray spectrum 
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Figure 3. Left: average OIR continuum during HS 1 . It is fitted with a power law mimicking the optically thin synchrotron from the compact jet = 
—0.63 ± 0.05), combined with a power law with a spectral index of ay = 2.73 ± 0.37. 

Right: ATCA+ISAAC+FORS2+XRT+PCA SED of GX 339-4 during HS1. The jet's optically thick power law (a sa 0.15, magenta dotted-line) is extended 
to the near-IR. The optically thin emission is roughly consistent, within the 3<r uncertainties (black dotted-lines) with the X-ray spectrum beyond about 
6 X 10 17 Hz, i.e. 2.5 keV. The accretion disc, derived from the X-ray fitting, is also superimposed. 



(Figure [4] right panel). To perform a more quantitative analysis, 
we used DISKIR combined with a jet near-IR power law to de- 
scribe the OIR/X-ray SED, but the best-fit, although statistically 
good, is rather unphysical. Indeed, the fraction of reprocessed X- 
rays / ou t is found to be between 0.1 and 0.5, and the ratio L c /Ld 
between the corona's and disc's luminosities as well as the temper- 
ature of the disc are unconstrained. This either means that the use 
of DISKIR is not relevant, or that disc irradiation combined with 
the optically thin emission of the jet is not sufficient to describe 
the OIR continuum of GX 339-4 during HS 1 . 

3.4 The average OIR continuum during HS2 

Figure [5] which displays the OIR spectra during HS1 and HS2, 
shows that the two spectra are very different. First, the HS2 spec- 
trum is clearly flatter in the near-IR, with a continuum lower in K 
but higher in H and J. Then, the flux drops dramatically and it be- 
comes fainter again from /. We can exclude any calibration issue 
as the explanation since the level reached in the optical matches 
very well to that of the J band. Moreover, we are not aware of 
any FORS2 and/or ISAAC instrumental problem that could have 
occurred, so that we are forced to conclude that this drop is real, 
although the phenomenon responsible for such peculiar continuum 
is difficult to understand. That said, we tested the straightforward 
explanation, i.e. the detection of the jet's spectral break, by fitting 
the spectrum with a model combining a broken and a simple power 
laws. We obtained a good fit, but the inferred slope (a $C — 2 be- 
yond the break at about 2.5 x 10 14 Hz) is steeper than expected 
for optically thin synchrotron (see Figure [6). This either implies 
a more complicated emission mechanism within the jet scenario 
and/or another origin for the OIR continuum during HS2. 



4 DISCUSSION 

In the SS, the OIR continuum can be best described by irradia- 
tion of the outer parts of the accretion disc. However, the fit to the 
data with an irradiated disc model implies a high level of repro- 
cessing, and an excess appears at lower frequencies. In the HS, not 
only do the near-IR jet-dominated spectra exhibit rapid variations. 



but at least one other process besides disc irradiation and optically 
thin synchrotron must be taken into account to explain the spectral 
continuum. 

4.1 GX 339-4 in the SS 

4.1.1 Enhanced irradiation of the accretion disc 

Irradiation of the accretion disc has long been thought to domi- 
nate the UV/optical emission of microquasars in the SS (see e.g. 
IVrtilek et al.lll990l : Ivan Paradiis & McClintocklll994h . and this is 
certainly the case for GX 339—4. However, the derived / out = 
0.05^0 02, although poorly constrained due to deredenning un- 
certainties, is at least twice the v alue previously measured for 
GRS 1915+105 dRahoui et alj201fjh . It would even be larger if we 
considered a higher value for the extinction along the LOS, (e.g. 
A v ~ 3.7. IZdziarski et al.lll998T) . We therefore believe that disc 
irradiation in GX 339—4 is enhanced, and two scenarios may be 
invoked: 

(i) IPringld Jl996b and lMalonev et alj d 19961) showed how non- 
axisymmetric radiation forces could induce warping of a flat disc. 
In such a configuration, the disc's surface illuminated by the cen- 
tral source would be larger, and the higher the degree of warping, 
the higher the illumination. It is therefore possible that the high 
level of X-ray reprocessing observed in GX 339—4 is the con- 
sequence, at least partially, of a ver y warped accretio n disc irra- 
diated by the inner regions (see e.g. lEsin et al]l2000l concerning 
GS 1124-68 and A0620-00). 

(ii) GX 339— 4's accretion disc is known to launch winds at 
large radii jMiller et al.ll2004h . Part of the X-ray and/or UV emis- 
sion could therefore be up-scattered by these winds back towards 
the disc to heat it up locally. This would also enhance the self- 
irradiati on of the disc, and would be added to the effect of warping 
(see e.g. lDubus etafl 19991) . 

4.1.2 The presence of a near-IR excess 

GX 339— 4's emission exhibits a near-IR excess with respect to 
the irradiated accretion disc emission. Before proceeding further 
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Figure 4. Same as Figure[3]but here, the spectral index of the optical power law, a y , was fixed to 2 to mimic the Raleigh- Jeans tail of the irradiated accretion 
disc. In this case, the optically thin emission (am = —0.74 ± 0.04) is not consistent, within the 3a uncertainties (black dotted-lines) with the X-ray 
spectrum. 



with our interpretation, we acknowledge that it could be due to dis- 
crepancies between the unknown spectral continuum of the stan- 
dard star HIP 68124 and the B2/3V synthetic model we made 
use of in place to flux-calibrate our data. That said, the excess's 
magnitudes in the 2MASS J, H, and Ks filters are 16.93 ± 0.34, 
15.70±0.22, and 14.76±0.15, respectively, approximatively sim- 
ilar to the lowest pub lished near-IR magnitudes to date (~ 17 , 
~ 16, and 15.2 ± 0.3. IChatv et al.ll2002l ; lBuxton & Vennesll2003l) . 
Considering that the source is in th e SS, when the radio activ- 
ity is quenched jFender et al.l 1 1999h . it is reasonable to assume 
that a compact jet contribution is at best marginal. We rather 
think that the near-IR excess is due to the secondary and, al- 
though a very simple modelling, we believe that the use of a black 
body is appropriate. The companion star has never been identi- 
fied, but based on the orbital parameters and lowest magnitudes 
in the op tical bands, it is thought to be a G, K , or M giant or 
subgiant dShahbaz et alj|200ll: [Hvnes et al.lr2004l : IZdziarski et al.l 
l2004l : lMunoz-Darias et al.l2008h . The temperature T B b ~ 3348 K 
that we derive i s consistent with a M star . Moreover, a mass ra- 
tio q ~ 0.125 flVlunoz-Dar ias et al.l l2008h and a semi-amplitude 
K, ~ 317 km s" 1 favnes et alj|2004l) lead to R,/a ~ 0.3 
and a sin (i) ~ 11.9 Rr* dZdziarski et all |2004|) . hence R* ~ 
3.6/ sin(i). The Rbb ~ 7.1 Rq v alue (for a minimum distance 
-Dbb ~ 6 kpc, Irlvnes et al.ll2003l) . not unusual for an M sub- 
giant, is therefore consistent if the inclination i ~ 30°, in agree- 
ment with previous estimates. This implies a ~ 23.8 Rq and 
Rout ~ 1.2 Rq, which again ar e consistent with the orbital pa- 
rameters of the system. However, IShahbaz et alj J200 lh measured 
r = 20.1 ± 0.1 in the off-state, and claimed that the undetected 
companion star was responsible for 30% of the flux at maximum, 
which gives r ~ 21.3. The inferred magnitude of the black body 
is r ~ 20.4. While this seems to invalidate the companion star 
hypothesis, such difference could actually be explained by orbita l 
variability (e.g. 0.4 magnitude in the / band. ICowlev et alj,l2002l ). 
companion star irradiation by the accretion disc, and, more im- 
portantly, by the discrepancies in the optical between a cool star 
spectrum and a black body. We therefore believe that this excess is 
indeed due to a M subgiant companion star, but further near- and 
mid-IR spectroscopic observations of the source in the off-state are 
necessary for confirmation. 



4.1.3 The HS1 case 

There is a strong contribution of the optically thin synchrotron 
emission of the jet to the OIR continuum. While this emission 
seems to extend to the X-ray domain, whether or not it means that 
the jet is responsible for the hard X-ray emissio n is discussed be- 
low. T his result is not new and was first given in lCorbel & Fender] 
Nevertheless, the use of spectroscopic OIR data here in- 
stead of photometric ones brings better constraints on the con- 
tinuum, especially in the optical. As already mentioned, accre- 
tion disc irradiation as modelled with the DISKIR model cannot 
alone explain the level of optical emission, as the derived fraction 
of reprocessed X-ra ys is 10-50%, against 5% in the SS. Indeed, 
Ide Jong etai] J 19961) found that a maximum of a few percent of the 
incoming photons were absorbed by the accretion disc, and it is 
hard to believe that disc warping and/or wind up-scattering could 
increase the fraction to such a level. Moreover, the spectral index 
of the optical power law may be too high compared to that ex- 
pected for the Raleigh-Jeans tail of a thermal process, although our 
data do not allow us to be definitive. We therefore believe that be- 
sides irradiation, there must be some other physical, perhaps non- 
thermal, process to account for the OIR excess in the emission of 
GX 339—4. A self-consistent model including a more complex jet 
and/or accretion physics is beyond the scope of this paper. Here, 
we mention two processes that we believe to be the best candidates, 
at least in the case of GX 339-4: 

(i) the first is pre-shock synchrot ron as described in th e frame- 
work of the jet model presented in iMarkoff et alj J2005b . It cor- 
responds to the emission of a magnetised reservoir of electrons 
coming from the accretion flow (similar to a corona), that remain 
in a thermal distribution at the base of the jet before they are ac- 
celerated into power laws. This region can produce a significant 
UV/OIR emission, whose intensity and peak frequency roughly 
depend on the region's size and the distance within the jet from 
where the acceleration starts, 

(ii) the second, recently proposed in IVeledinaetall i l201ll) . also 
involves synchrotron radiation from electrons located in a corona 
close to the BH. Nonetheless, the population has a non-thermal 
distribution. This region produces a strong OIR emission, whose 
intensity and peak frequency mainly depend on its size, and an X- 
ray/OIR anti-correlation is found. Combined with the time-lagging 
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Figure 5. Comparison between the GX 339— 4's OIR deredenned spectra, 
obtained during HS1 (gray) and HS2 (black). 



X-ray/OIR correlation due to X-ray reprocessing within the disc, 
this model can explain the pattern obse rved in the optical/X -ray 
cross-correlation function of GX 339—4 dGandhi et alj|2008h . 



4.1.4 The HS2 case 

4.2 GX 339-4 in the HS 

We are unsure of the phenomenon responsible for such peculiar 
continuum. A possible interpretation is that it is due to the same 
process that might be responsible, in addition to disc irradiation, 
for the optical excess during HS1. Indeed, we may be witness- 
ing a shift of this component to the lower frequencies, where it 
peaks in H or J. This would explain the flatter slope of the opti- 
cal excess, more consistent with disc irradiation. Such shift is pre- 
dicted for the two p rocesses described in Sect. 4.2.1 . For instance, 
iMaitra et alJd20091) , using the jet model presented in lMarkoff et al.l 



showed how changes in the emitting region's size and/or 
the distance from which acceleration starts led to pre-shock syn- 
chrotron peaking either in the UV/optical or the near-IR (Figure 4 
in their paper). However, another explanation could be a more 
c omplex compact jet's sp ectral continuum, such as those presented 
in lPe'er & Casellal | |2009|) . In their paper, the authors take adiabatic 
losses into account, both for power law and Maxwellian distribu- 
tions. Depending on the value of the magnetic field, they show 
that this gives a large sample of continuum shapes. In particular, 
Figures 6, 7, and 8 in their paper, which displays profiles corre- 
sponding to a Maxwellian distribution for different magnetic field 
values, could correspond to the compact jet emission during HS2. 

4.2.1 Understanding the rapid near-IR continuum variability 

During HS1, the variability must be mainly driven by rapid fluctu- 
ations in the compact jet, and it is interesting that this phenomenon 
has almost no effect on the slope. This is a hint that the optically 
thin synchrotron emission is rapidly moving upwards or down- 
wards. This may imply rapid changes in the l ocation - in flux an d 
frequency - of the undetected spectral break dGandhi et alj|201lh , 
i.e. significant variations at timescales at least as low as 20 s of 
the physical conditions at the base of the jet, such as the magnetic 
field and/or the base radius. Moreover, the fact that during HS1, 
the optical, J, H, and K spectra, although non-simultaneous, are all 



Figure 6. Average OIR continuum during HS2. It is fitted with a broken 
power law («ir,i ~ —0.47 and ctiR2 ~ —2.94), combined with a power 
law with a spectral index fixed to ay = 1.5. Clearly, aim is too steep for 
optically thin synchrotron radiation. 



consistent with each other and that the extension of the near-IR 
jet power law to the high energy domain roughly describes the 
3-50 keV RXTE spectrum is surprising. It means that the opti- 
cal and near-IR sub-second flickering observed in the HS (se e.g. 
iGandhi et al ]|20ld : ICasella et alfeOld) are smoothed out and that 
we observe the average behaviour of the compact jet beyond about 
a few hundred seconds. This is in good agreement with the exis- 
tence of tight near-IR/X-ray correlations in microquasars in gen- 
eral and GX 339 —4 in particular (see e.g. iGleissner et alj|2004l 
ICoriat et"al]|2009l). In any case, in a recent study of Cygnus X— 1 
l Rahoui et al. 201 if) , we argue that the extension of the optically 
thin synchrotron emission from the spectral break we detect in the 
mid-IR to the 7— ray domain matches pr etty well to the ten-y ear 
averaged polarised emission mentioned in lLaurent et alj J201 ll) . In 
contrast, the contribution of the compact jets to the X-ray spec- 
trum (^Xr£/PCA+HEXTE) is found marginal, a hint that it prob- 
ably comes from a hot corona. In this context, a confirmation of 
the domination of compact jets in the GX 339— 4's X-ray emission 
could mean that both microquasars belong to two distinct families, 
perhaps differentiated by the mass-accretion rate level and/or lu- 
mino sity dRussell et al.l2010i : ICoriat et alj201 U lSobolewska et all 
120111) . 

Considering that we do not clearly understand the origin of 
the near-IR emission of GX 339—4 during HS2, it is difficult to 
speculate about the observed variations. The continuum is more 
variable in H and J than during HS1 and, although K exhibits a 
higher maximum-to-minimum flux ratio, it only stems from three 
sub-spectra and / max //mm = 1.15 ± 0.04 otherwise. This be- 
haviour could mean that the variations are not driven by the same 
processes, or that the compact jet contributes more to the H and J 
bands, in other terms that the spectral break is located in the near- 
IR. It is also puzzling that the discrepancies observed in the OIR 
between both HS have no equivalent in the RXTEfPCA spectra. 
The reason could be that the X-ray data were obtained more than 
half a day after the OIR ones, i.e. an event could have occurred 
while GX 339—4 was being observed in the OIR before returning 
to the "normal" during the RXTE observations. But it could also 
mean that the hard X-ray emission is due to another component, 
and whatever happens to the jets has no consequence in the high 
energy domain, at least not at these timescales. 
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5 CONCLUSION 

We presented a comprehensive multi-wavelength study of 
GX 339—4 in the SS and HS, with a focus on the OIR spectral con- 
tinuum. We showed that the SS is dominated by enhanced irradia- 
tion of the accretion disc, with an additional contribution from a M 
subgiant companion star. In one of the HS, the OIR emission stems 
from the compact jet, disc irradiation, and maybe synchrotron from 
the inner regions. The second HS spectrum exhibits a peculiar con- 
tinuum which could be a hint of a more complex emission. In both 
cases, the near-IR spectrum is strongly variable on timescales at 
least as short as 20 s, which we suggested is a manifestation of 
rapid variations of the magnetic field and/or nozzle radius. Our 
results highlight the complexity of the OIR emission of a micro- 
quasar like GX 339—4, and a good understanding of the source's 
behaviour in this spectral domain requires further in-depth spectro- 
scopic observations as well as theoretical modelling. In particular, 
the use of an instrument allowing simultaneous integrations in the 
whole OIR range, such as ESO/X-shooter, would be of great help 
to investigate the variations - at long and short timescales - of the 
continuum in a more robust way. Finally, this strengthens the need 
of the lames Webb Space Telescope (JWST) for microquasar stud- 
ies, as the unmatched sensitivity of NIRSpec and MIRI will allow 
a simultaneous spectroscopic coverage from 1 to 30 fim. 
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Figure Al. Synthetic spectra used for the near-IR spectroscopic standard 
stars in this study. Their measured B, V, J, H, Ks flux densities are super- 
imposed. 
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Figure A2. Synthetic spectra used for the optical spectroscopic standard 
stars in this study. Their measured B, V, J, H, Ks flux densities are super- 
imposed. 



APPENDIX B: COMPARISON OF GX 339 4'S OPTICAL 
SPECTRA DERIVED FROM THE TWO METHODS 

Figure IbTI displays a comparison, for each observation, between 
the optical spectra of GX 339—4 obtained using either the ESO 
spectro-photometric calibration or Kurucz templates (left panel), 
as well as a zoom-in the optical/7 overlapping region (right panel). 
The use of synthetic spectra is necessary to extend the wavelength 
coverage and cover the gap with the near-IR. 



APPENDIX A: FLUX-CALIBRATED SYNTHETIC 
TELLURIC SPECTRA 

In the near-IR, the lack of spectro-photometric standard stars im- 
plies the use of synthetic spectra for flux-calibration purpose. Fig- 
ure |Al] displays the Kurucz spectra we used in the near-IR, and 
Figure lA2l displays the ones we used in the optical in order to ex- 
tend the wavelength coverage beyond 1.02 /im. 
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Figure Bl. Observed optical and J band spectra of GX 339—4 during Obs. 1 (top), 2 (middle), and 3 (bottom), full-scale (left) and zoom-in on the overlapping 
region (right). The optical spectra derived using a synthetic standard star template are displayed in red and the ones using the ESO spectro-photometric 
calibration in green. The J band spectra are displayed in blue. Note that the telluric absorption features are still present in the green spectra. The red spectra 
displayed here are rebinned to a third of their real resolution to improve the visibility of the optical/7 overlapping region. However, only the original ones 
were used for SED fitting. 
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